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Abstract 



Aims. We use a specific form of the interpolation function in the MOND formalism, which optimally accounts for the 
internal structure of dwarf spheroidal (dSph) galaxies, to explore the consequences it has on the scaling relations seen 
in these systems. 

Methods. The particular form of the interpolation function we used leads to a law of gravity that does not degrade the 
good fit of the MOND proposal on galactic scales, and in fact, slightly improves the accordance with observations at 
dSph scales. This formalism yields a good description of gravitational phenomena without the need of invoking any still 
'^'i undetected and hypothetically dominant dark matter, in the weak field regime probed by local dSph galaxies, 

'^s ' Results. Isothermal equilibrium density profiles then yield projected surface density profiles for the local dSph galaxies 

I in very good agreement with observational determinations, for values of the relevant parameters as inferred from recent 

• , observations of these Galactic satellites. The observed scaling relations for these systems are also naturally accounted 

i-J^ ■ for within the proposed scheme, including a previously unrecognised correlation of the inferred mass-to-light ratios of 

local dSph's with the ages of their stellar populations, which is natural in modified gravity schemes in the absence of 
JL I dark matter. 

t*^ , Conclusions. The results shed some light on the form that the MOND interpolating function may have in the most 

■ challenging regime, which occurs at moderate accelerations and intermediate mass-weighted lengths. 
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> 1. Introduction et al. (2007), Famaey et al. (2007), Gentile et al. 2007), 

^ . . Tiret et al. (2007) and Sanchez-Salcedo (2008). 

. If one wants to explore options where the existence of 

' dark matter is not a necessity, one must consider vari- „, r i i- i • i i i ^ ^ t 

... • ,1 1 r 1 i- J i 1 i- 1 he range ot galactic dynamical problems treated With 

ations m the law ol gravity on galactic and extra-galactic h^/^att^ r i- i j.jjr J.^ n i. j 

, 1 J, °. 11 J, , r the MOND formalism has extended Irom the hrst-order 

■«H- . scales where dark matter has been proposed to account tor i rr ^ r ^ i- i i t 

, . , , , J . rrii 1 , I J ■ 1 1 1 ■ gravitational enects ot rotation curves and velocity disper- 

' the observed dynamics. Ihe best-studied such proposal is . ^ ^ . , r V^i 

• njrr^j-c J AT J. • T-v • ^n/r/^ATT-.\ 1 A. • c sion measurements to cover a wide range ot more subtle 
the MOdihed Newtonian Dynamics (MOND hypothesis ol ttt-^i • i • t -i- 5 ■ r ^ i 

O n/r-1 /1noo^ 1-11 \ u j- j- r problcms. With their role m limiting the sizes ol satel- 

Milgrom (1983), which has been shown to account tor the ^ ^ ^ ^■ t ■ ■ h^^att^ i -i- i- ^ i 

^ r-11-/ oj pn/rri I. lites and establishing m MOND escape velocities tor satel- 

> rotation curves ot spiral galaxies (e.g. banders & McGaugh , . ^ ■ j_ ^ ^ ^ ^ ^- n ^ ^ 

onnoN J- 11 ■ i- -Li rrTii tt i w iites Or galaxics subiect to an external acceleration held, 

• ^ , 2002), naturally incorporating the iuUy-i'isher relation, see i r ^ j- j o i o i j c tt j 

r ^ ^„■^ fnnno f J r • tidal torccs were studied by Sanchez-Salcedo & Hernandez 

■ tor example Milgrom (2008a) and relerences therein. /n^^^^-T^ ttt ^ i /nnnn\ ^ ^ ^ ^ ^^ i r 

: Al 1 11 (2007). Wu et al. (2008) calculated the escape velocity lor 

^ Abundant recent publications of the velocity disper- ^^jj^y ^^^^^ ^^^^ ^^q^j3 ^^^j^ ^^^^^^^ 

. . > sion measurements of stars m the local dwarf spheroidal conclude that the LMC appears as bound from both points 
(dSph) galaxies, the extended and flat rotation curves of ^^^^^ ^pi^^ recently determined high proper 
spiral galaxies, the high dispersion velocities of galaxies m ^^^^j^^^ ^^^^ ^^-^^^^ ^^^^ ^ligl^^ly ^^^^^ 
clusters, the gravitational lensmg due to massive clusters of ^q^^ prescription. Sanchez-Salcedo et al. (2008) looked 
galaxies, and even the cosmologically inferred matter con- ^j^^ thickness of the extended HI disk of the Milky Way 
IT^J?^ universe, are all successfully modelled using ^^^^ ^^^^^^ ^ somewhat better fit to ob- 
MOND, not as indirect evidence of a dominant dark mat- ggrvations in the MOND theory. Sanchez-Salcedo et. al 
ter component, but as direct evidence of the failure of the (2OO6) and Nipoti et al. (2008) examined the problem of 
current Newtonian and general relativistic theories of grav- dynamical friction in dSph galaxies, comparatively assum- 
ity, m the large-scale or low-acceleration regimes relevant to Newtonian gravity or MOND, with decay timescales for 
the above. Some recent examples of the above are Milgrom globular clusters being somewhat shorter in MOND. This 
& Sanders (2003), Sanders & Noordermeer (2007), Nipoti ^ potential problem, unless one allows large initial orbital 
radii for the observed globular clusters beyond the current 
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Going to cosmological scales, Skordis et al. (2006) have 
studied the cosmic microwave background, Halle et al. 
(2008) looked at the problem of the cosmic growth of struc- 
ture, Zhao et al. (2006) studied gravitational Icnsing of 
galaxies, and Angus et al. (2007) and Milgrom & Sanders 
(2008) studied dynamics of clusters of galaxies. All of the 
above find the MOND description of the problem to be a 
viable option within the observational errors of the relev- 
ant determinations, despite the need to include some unseen 
mass on galaxy cluster scales. Options for the above include 
the cluster baryonic dark matter proposed by Milgrom 
(2008b), or the inclusion of neutrinos, whether active as 
studied by Sanders (2003) and Angus et al. (2007) or sterile 
as proposed by Angus et al. (2010), always within plaus- 
ible expectations. Also, in spiral galaxies, although within 
observational errors, rotation curves have been reported by 
Gentile ct al. (2004) and CorbclH & Salucci (2007) to be 
somewhat less well-fitted assuming the MOND prescription 
than when using cored dark matter haloes. 

Observationally, it appears that MOND fares as well 
as dark matter in accounting for measured dynamics, in 
all but the smallest scales, the case of local dSph galax- 
ies remains the most controversial. Published studies using 
MOND, which typically assume the deep MOND regime, 
sometimes find a different value for the acceleration scale of 
the theory for different dSph galaxies, when fitting detailed 
dynamical models to the data (Lokas 2001), or that M/L 
ratios remain higher than those of stellar populations, mak- 
ing dark matter necessary even under the MOND formula- 
tion, e.g., Sanchez-Salcedo & Hernandez (2007). Still, this 
last point remains somewhat controversial. Angus (2008) 
suggests that the MOND proposal will hold for local dSph 
galaxies, once the effects of tidal disruption are adequately 
incorporated. 

Regarding corresponding theoretical developments, nu- 
merous alternative theories of gravity have recently ap- 
peared (Bekenstein 2004, Sanders 2005, Sobouti 2007, 
Bruneton & Esposito-Farese 2007, Zhao 2007, Arbey 2008), 
now mostly grounded on geometrical extensions of gen- 
eral relativity and field theory, which lead to laws of grav- 
ity in the Newtonian limit that in the large-scale or low- 
acceleration regime reduce to the MOND prescription fit- 
ting formula. 

The MOND theory is characterised by (i) a low accel- 
eration regime where dynamics mimic the presence of dark 
matter, (ii) a high acceleration regime where Newtonian 
gravity is recovered, and (iii) a somewhat ill-defined trans- 
ition region. Variants of MOND are defined by the form of 
the interpolation function, /^(.t), which mediates the trans- 
ition between the two limit regimes, sec Eq.(l). 

We show how the MOND prescription with a suitable 
interpolation function, such as the one used by Bekenstein 
(2004) and used by Famaey & Binney (2005), can be writ- 
ten as the addition of both low acceleration limit MOND 
and Newtonian contributions to the acceleration on all 
scales. The large differences in scales and magnitudes of the 
two acceleration terms ensure that the addition does not 
spoil the good match with observations on galactic scales. 
We are trying to look for a suitable MOND interpolation 
function for the local dSph galaxies, but not necessarily one 
that works on all scales. 

This, as it is the regime of the local dSph galaxies, the 
one which has proven most difficult for MOND, and the 
one where the high-quality recent observations now avail- 



able allow for a relatively clean test. We find that this in- 
terpolation function not only yields acceptable M/L values 
for all dSph galaxies without the need for any dark matter 
in all cases, h\it also provides a natural explanation for all 
the scalings seen in these Galactic satellites. An interesting 
correlation between the M/L ratios inferred and the ages of 
the stellar populations in local dSph galaxies is found, and 
naturally accounted for in modified gravity schemes where 
stars alone account for the gravitational force. 

However, on solar system scales, Sanders (2006) has 
shown this particular interpolation function to be incom- 
patible with limits on the variations of Kepler's constant. 
Thus, the good accordance we find with observations for 
dSph's must be understood as evidence in favour of the 
particular interpolation function we test, at the weak field 
limit relevant for these systems. What a full interpolation 
function for MOND should be is a more complicated issue. 
One, however, which has to be addressed bearing in mind 
the results presented here. 

Section (2) gives a brief summary of the interpolation 
function first explored by Bekenstein (2004) and gives equi- 
libriiim isothermal configurations. These are then used in 
section (3) to model the local dSph galaxies and obtain 
M/L ratios, which are found to be consistent with those 
typical of stellar populations, at the ages of the different 
dSph systems. Section (4) includes an exploration of the 
most conspicuous scaling relations for local dSph galaxies, 
finding that Bckenstein's interpolation function naturally 
accounts for all features found. Our conclusions are sum- 
marised in section (5). 

2. A particular form of the interpolating function 
in the MOND prescription 

In terms of the acceleration g felt by a test particle, the 
MOND proposal is (see e.g. Milgrom 2002) 

= (1) 

where Qj^ is the acceleration assigned by Newtonian gravity, 
ao an acceleration scale of the theory, and the interpolat- 
ing function /i(a; := g/aa) is an unspecified function that 
reduces to unity for large values of its argument, and to 
its argument, for small values of it. In this way, Newtonian 
gravity is recovered for large values of g, and in the deep 
MOND regime one obtains 

g = (ao.gAf)^^^ • (2) 

Strictly speaking, MOND is a theory defined only at 
the limit values for the acceleration. The value of the con- 
stant flo has been reasonably determined by calibrating the 
deep MOND regime through observations of rotation velo- 
city curves of large spiral galaxies (e.g. Sanders & McGaugh 
2002). However, the details of the transition region between 
the two regimes (i.e. the function /x) have proven harder to 
establish reliably. This is in part natural, as it is not a scalar 
parameter, but a functional dependence that one is look- 
ing for here. A number of variants for this function have 
been proposed (e.g. Famaey & Binney 2005, Zhao 2007, 
Sanders & Noordermeer 2007), but apparently suitable dy- 
namical systems in the transition region arc hard to come 
by. Astrophysically, to first approximation one tends to find 
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either systems where no dark matter is needed (stellar and 
planetary systems, globular clusters, vertical dynamics of 
disk galaxies, galactic bulges, and elliptical galaxies), or 
systems where dark matter is massively dominant (dSph 
Galactic satellites, rotation curves of spiral galaxies, galaxy 
groups and clusters, and cosmological observations). Also, 
the details of any empirically inferred transition region, 
/x(a;), are sensitive to the details of the baryonic system one 
is looking at, such as the assumed mass-to-light ratios, gas 
dynamics, or orbital anisotropy of stars, most of which play 
a marginal role when in the deep MOND regime. As poin- 
ted out already, this ill-defined transition is cumbersome to 
handle in attempts to reproduce the MOND phenomeno- 
logy through simple GR extensions e.g. Sobouti (2007) or 
Capozziello et al. (2007). 

Bekenstein (2004) shows that his relativistic extension 
of TeVeS in the appropriate non-relativistic limits yields 
the interpolating function 



VI + 4x - 1 



(3) 



This particular interpolating function converges to the right 
limits as x — > 0, cxi and has a very peculiar property. 
Direct substitution of equation ^ into the absolute value 
of relation ^ yields 



9 ^ 9n + {a()9N ) 



(4) 



This equation can be thought of as a generalised gravity 
recipe described by the addition of two terms, the first 
a standard Newtonian acceleration term and the second 
the MOND limit acceleration term. Seen in this way, equa- 
tion Q changes its acceleration behaviour limiting cases 
to a scale limiting behaviour. Indeed, for the case of a test 
particle on a gravitational field produced by a central mass 
M, located at a distance distance R from it, equation Q 
can be written as 



9 



GM [Ga^f^ Ml/2 



(5) 



i?2 R 
where G is Newton's constant of gravity. Seen in this way, 
equation ([5]) converges to Newtonian gravity for sufficiently 
small i?/Mi/2's and reproduces the MOND strong limit for 
sufficiently large R/M^/'^^s, and so the acceleration lim- 
its are now "scale- weighted" limits. It is worth noting 
that equation ([S]) is certainly the simplest modification to 
Newtonian gravity once the weak acceleration MOND re- 
gime is known. 

Let us now rewrite equation Q in such a way that 
the scale weighting becomes clearer. To do so, recall that 
gw = — Vt/fAT, where (f>N represents the standard Newtonian 
potential and g = —V4>, where (j) is the scalar potential of 
the gravitational field. With this, equation (g]) can be writ- 
ten as 

9 = 9N {I + x) , (6) 

where x ■= (ao/G)!/^ {R/W^/^). In what follows we use the 
standard values of oq = 1 x lO^^cms^^ (Milgrom 2008a) 
and G = 4.5 X IO^^^A/q ^ s"^ kp^s^ -^^ ^^^j^g suitable for 
galactic applications. 

As already mentioned, we are interested in the form of 
Ijl{x) at scales of the local dSph galaxies. That the partic- 
ular pl(x) we are testing is incompatible with solar system 



dynamics implies that a complete /i(a;) MOND function is 
probably complex, that is, if one wants to treat the MOND 
formalism as more than just an empirical phenomenological 
description of gravitational physics on galactic scales. On 
the other hand, some researchers (see e.g. Bekenstein 2006 
and references therein) claim that the Pioneer anomaly 
can be explained by the MOND hypothesis. This assump- 
tion can prove wrong if uneven thermal radiation in the 
spacecrafts is found (Toth & Turyshev 2009), which could 
possibly come from the flyby anomaly (Turyshev 2009). In 
other words, the MOND formalism does not appear to be 
relevant on solar system scales, where General Relativity 
and Newtonian Gravity have proven correct (e.g. Turyshev 
& Toth 2009, Anderson et al. 2002). Still, as discussed by 
Milgrom (2009), it is possible that current tests on solar 
system scales cannot reach definitive conclusions on the 
MOND interpolating function. 

For globular clusters, with M = IO^-^Mq, i? = 2 - 10 x 
10~'^kpc, we get values for x of between 0.1 and 0.01, and 
the correction becomes smaller than the errors and uncer- 
tainties in the observational determinations for the values of 
radii and masses for globular clusters. For elliptical galaxies 
and bulges, with masses going from about 10^ to IQ^^Ad^ 
and radii of between 0.5 and 10 kpc, x is about 0.1. We 
thus see that the correction to gravitational dynamics due 
to the proposed inclusion of a second term in equation ([6]) is 
small enough to have remained undetected in galactic sys- 
tems where no gravitational anomaly is found and where 
dynamics are consistent with Newtonian gravity, in the ab- 
sence of any dark matter. 

On the other hand, for systems where the presence of 
dark matter is inferred, given that this is always required to 
be dominant, adding the Newtonian term onto the MOND 
proposal generally provides a negligible contribution. For 
example, for the Galactic disk at the solar radius, M = 
5 X IG^^Mq and R — 8.5 kpc give x = 2, consistent with 
an inference of about 50% dark matter within the solar 
circle. In going to the outskirts of the Milk Way, we go to 
R = lOOkpc, hence x — 20. The system is either totally 
dominated by dark matter, or it is in a regime where the 
second term in equation ([6]) almost fully determines the 
dynamics. We see that the distinct power-law dependences 
ensure that the Newtonian term completely dominates at 
small values of x, while the opposite holds for large values 
of X, with a necessarily narrow transition region. 

From the form of equation ([5]) , it would be tempting to 
add the following term of the ^ 1/R series, a further (per- 
haps positive) constant term, which would result in an x^ 
additive term in equation ([B]). If chosen suitably small, for 
the same reasons as given above, it would have no meas- 
urable effects on all but the largest scales, perhaps as a 
tool to model the cosmological constant. This extension 
certainly does not correspond to MOND, since no interpol- 
ation function fi{x) can be constructed. However, in such 
terms, equation © might be interpreted as a series expan- 
sion of a more fundamental underlying gravity law, only the 
first terms of which we have begun to appreciate empiric- 
ally as observations probe increasingly higher x regimes, 
typically corresponding to increasingly larger scales. The 
introduction of the constant term in equation ([5]) and its 
calibration from cosmology, however, must be done within 
a generalised GR framework. 
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2.1. Equilibrium configurations 



To test at the regime which has resulted most troublesome 
for the MOND formalism, we go to the now very well stud- 
ied dSph galaxies of the local group. To compare against 
local dSph galaxies, we require the derivation of equilibrium 
configurations for a population of self-gravitating stars. We 
begin by noting that the validity of Newton's theorems for 
spherically symmetric matter distributions, that the con- 
tribution to the force felt by an observer due to external 
shells vanishes, and that the contribution of all shells in- 
terior to the observer is equivalent to concentrating all 
mass interior to the observer at the centre, also holds under 
the present proposal. The two theorems depend on the nu- 
merator of the Newtonian term in equation ([S]) for a fixed 
solid angle fraction of a thin shell scaling with the second 
power of distance to the shell, as the denominator does. For 
the second term in equation ^ , the numerator scales with 
^y{R'^) = R, as does the denominator, assuring the validity 
of Newton's theorems for spherical mass distributions (e.g. 
see Bekenstein & Milgrom 1984). 

We now write the equation of hydrostatic equilibrium 
for a polytropic equation of state p = Kp'^: 



GM{r) (ao/G)'/'M(r)i/2 



(7) 



dr 

{4TTr'^)~^dAI{r) /dr, where M{r) is the mass of 



Since p 

the configuration at a radius 



conditions, K 
be written as 



and going to isothermal 



a with 7 = 1, the previous equation can 



\ dr I dr'^ r 



GM{r) {ao/Gf'^M{rf/'^ 



In the above second-order differential equation for 
M(r), the constant isotropic velocity dispersion for the pop- 
ulation of stars has been written as a. We have taken ad- 
vantage of the correspondence between hydrostatic equi- 
librium polytropic configurations and self gravitating stel- 
lar systems, which does not depend on the particular form 
of the Poisson equation, only on the validity of the iso- 
thermal condition for the stellar population (see e.g. Binney 
& Tremaine 2008). Taking initial conditions M{r) — >■ and 
dM/dr = Anr'^po for r — )■ 0, a constant central density po, 
we can now solve equation (8) through a numerical finite 
differences scheme, subject only to two input conditions, a 
value for a and a value for po- 

Notice that neglecting the first term on the right hand 
side of equation (8), the limit behaviour for M{r) for large 
values of r when the weak limit of the MOND formal- 
ism dominates, is characterised by a finite total mass. In 
contrast to what happens in Newtonian systems where 
isothermal configurations have an infinite total mass, iso- 
thermal self-gravitating configurations will be naturally 
bound in mass, under the proposed gravity law. As a result, 
there will also be a well-defined and finite half-mass radius, 
^^hm, to characterise the resulting equilibrium configura- 
tions. Also, as is the case using the MOND prescription, 
the total mass of the configuration is expected to scale with 
a . From the dimensional mass scale < M >~ (ao/G) \ 
one is led to expect an analogue to the TuUy- Fisher relation 
for spheroidal galactic systems, down to the dSph regime. 



-1.5 -1 



log(R) 



Figure 1. A sample isothermal equilibrium density pro- 
file for po = 0.1Mqpc~^ and a = 7km/s, resulting in a 
total mass of 3 x IO^Mq and a volume half-mass radius of 
0.39kpc. The density is measured in units of Mqpc~^ and 
the distance in units of kpc. 

Figure (1) gives a plot of log(p/M0 pc"'^) vs log(i?/kpc) 
for a sample numerical solution to equation (8) for a = 
7kms~^ and po = O.IMqpc"'^. We obtain a well-defined 
resulting total mass, Mtot, of 3 x IO^Mq, and a final volume 
half-mass radius of 0.39 kpc. We see the expected decay of 
the density profile at large radii, slightly faster than r~^. 
It is this type of density profiles that is used in the next 
section to model the local dSph galaxies. 

3. Local dwarf spheroidal galaxies 

As discussed above, once values for the velocity dispersion 
and the central density are given, numerically solving equa- 
tion (8) yields the full equilibrium density profile. To now 
model dSph galaxies, we take values for the reported ve- 
locity dispersions for these systems, which thus fix one of 
the two parameters of the model. A resulting density pro- 
file is then integrated along one dimension to yield a surface 
density profile, from which the projected half-mass radius is 
measured. This provides a second constraint; we then vary 
the value of the input central density to ensure that the 
resulting half-mass radius i?hm matches the reported half- 
light radius i?hi of a given dSph. Observational determin- 
ations of the velocity dispersion and the half-light radius 
of a certain dSph galaxy then fully determine the model. 
Dividing the total mass for the resulting model Mtot by the 
reported total luminosity Ltot of a dSph, yields a mass to 
light ratio for the model. Values of the half-light radii and 
total luminosities were taken from Gilmore et al. (2007), 
who summarise published results from the references given 
in Table (1), except for Ursa Minor, where we took the 
revised value for Ltot from Palma et al. (2003). 

Although other stellar systems show velocity dispersion 
profiles that generally decay as the radial coordinate in- 
creases, e.g. globular clusters, the case for observed dSph 
galaxies is different, with these systems showing essentially 
fiat velocity dispersion profiles, e.g. Gilmore et al. (2007). 
In cases where some drop in the velocity dispersion profile 
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Galaxy 


(J |kni / s] 


fknrl 


Lt„t X lO^Lr^ 

-'--'tot -'-'Cl) 


(MIL) 4 


(M/L) 


< V > 
^ A ^ 


Age of Youngest 
Component [Gyr] 


Carina 


7± 1.8 


0.290 


4.3 


O.D„2.9 




7 


3 


Draco 


8± 1.5 


0.230 


2.6 


AO q+29 


170+13.9 

^^'•"-8.9 


4.1 


10 


Leol 


8 ± 1.2 


0.330 


48.0 


n -70.65 
' -0.3 


r.U_o.44 


5.7 


2 


Sextans 


7± 1.0 


0.630 


5.0 


q r,+5.3 


6.3t*:^ 


13.4 


(2-6) 


Fornax 


12 ± 1.3 


0.400 


150.0 


1 4+0.45 
J^-*-0.35 


i.4ir4 


> 3.4 


(2-3) 


Sculptor 


9.5 ± 1.7 


0.160 


22.0 


Q 7+2.2 
■3. ' _i.4 


4+2.1 


> 2.2 


> 5 


LeoII 


6 ± 1.4 


0.185 


7.0 




0+2.4 


5.5 


6.5 


Ursa Minor 


8 ± 2 


0.300 


5.8 


r + 6.5 
0-°-3.6 


8.01'-? 


> 5.3 


12 



[M/L) A gives the values for the mass to light ratios calculated by Angus (2008), and M/L those in this study under the 
proposed model. Total luminosities (in the V band) and half-light radii are from Wilkinson et al. (2006), Wilkinson et al. (2004), 
Koch (2007a), Kleyna et al. (2004), Walker et al (2006), Mateo (1998), Coleman et al. (2007) and Irwin & Hatzidimitriou (1995), 
for the galaxies in the Table, in the order given, as summarised in Gilmore et al. (2007). Velocity dispersions are from adjusting a 
constant value to the data of Angus (2008). < X > gives the average value of the parameter x defined in section (2), which 
measures the relative relevance of the Newtonian and the low acceleration limit of the MOND formalism terms in equation (6). 
The final column gives an estimate of the age of the youngest stellar population present in each of the systems, for Carina, Ursa 
Minor, Leol, and LeoII from Hernandez et al. (2000), for Fornax from Coleman & de Jong (2008), for Draco from Aparicio et al. 
(2001), for Sextans from Lee et al. (2003), and for Sculptor from Babusiaux et al. (2005). 



is observed, this typically occurs towards the edge of the 
galaxy, affecting only a very small percentage of the total 
mass of the system. The constancy of these observed ve- 
locity dispersion profiles validates the use of equation (8) 
under isothermal conditions for modelling local dSph's. We 
have not at this point attempted a more detailed model- 
ling considering non-isothermal conditions, for example in- 
cluding orbital anisotropy, although such models are now 
common in models of dSph galaxies under Newtonian (e.g. 
Lokas 2002) or MOND frameworks (e.g. Angus 2008) and 
could be included in subsequent analysis. We point to the 
recent work by Gilmore (2007), in which isothermal con- 
ditions are used to model local dSph galaxies assuming 
Newtonian gravity. 

Values of a were taken from adjusting a constant level 
to the a profiles of Angus (2008), which are flat to a very 
good approximation, with the exception of Ursa Minor 
that shows a significant decrease with increasing radial dis- 
tances. This galaxy also shows internal structure in phase 
space at small radii, which might partly account for the 
steep increase in a seen by Angus (2008) towards the centre. 
We have taken a value for a representative of the situation 
at around the half-light radius. If a strong radial depend- 
ence of a in Ursa Minor were confirmed, however, this sys- 
tem would have to be excluded from the present sample, as 
it would then conflict with the simple isothermal modelling 
we are performing here. 

Table (1) gives our results for the 8 best-studied local 
dSph galaxies. The values we obtain for M/L are mostly 
comparable to those obtained by Angus (2008) for the 
same systems, but a slight systematic decrease is evident, 
consistent with having included a further force term, the 
Newtonian component, which reduced inferred AL/ L ratios 



even further. In general, the values of the parameter < x > 
we found, calculated from the resulting model total mass 
as {aQ/GYl'^{2^l'^R^y/Mll^), explain why the addition of 
the Newtonian term is required to fully account for the dy- 
namics. Further, the wide range of values of < x > shown 
in the table explains why a different ao has sometimes been 
found for various dSph's, when fitting dynamics using only 
the weak acceleration limit of the MOND formalism, e.g. 
Lokas (2001). 

The case where the effect of the second term in equa- 
tion (|6]) is greatest is that of Draco, where the compact con- 
figuration and relatively high a values allow the Newtonian 
term to figure somewhat more, leading to a strong decrease 
in the inferred M/L of between 8.1 and 30.9, compared to 
the Angus (2008) results of M/L between 22.6 and 72.9. 
With the possible exception of Fornax and Sculptor, cases 
where the reported King core radius serves only as a lower 
limit for the projected i?hi, the value for the parameter 
< X > for Draco is the lowest in the sample. This naturally 
explains why it is here that our results for M/L differ most 
and show a significant reduction when contrasting with the 
results of Angus (2008), who takes a fi function that rap- 
idly converges to the low acceleration limit of the MOND 
prescription, as this is the parameter that determines the 
relative importance of the Newtonian and low acceleration 
limit MOND terms in equation 

The 1(7 confidence intervals of our results are only lower 
bounds, as their calculation includes only uncertainties in 
the adopted values of a and not observational errors in total 
luminosities and half-mass radii. Also, freedom in the light 
profile functional fitting would increase our confidence in- 
tervals slightly, e.g., the reported i?hi values for Ursa Minor, 
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Figure 2. Comparison of our projected surface density pro- 
file for an cqiiilibrium isothermal solution to equation (8), 
having input a and projected i?hm as the observationally 
determined values for LeoII, dashed curve. For comparison, 
we also give the star counts surface density profile for Leo 
II of Coleman (2007). Both profiles have been normalised 
to the same total luminosity. 



Fornax and Sculptor are only lower bounds, because the 
quantity available for those galaxies is the core radius from 
King profile fits. Still, this last is only a second-order effect, 
since our inferred M/L values are not very sensitive to the 
values of Rhi used, provided this does not change by a large 
factor (see the following section). 

As expected, the main determinant of the resulting total 
mass of a numerical solution of equation (8), is the input 
a. We find, as expected from the situation of the low ac- 
celeration limit of the MOND prescription, a strong correl- 
ation, with the resulting scaling for all galaxies being well 
described by 



2\2 



Mtot = 



(3.36a') 
aoG 



(9) 



This could be seen as an extension of the Tully-Fisher re- 
lation down to the smallest galactic scales, a natural con- 
sequence of the model being explored. 

Considering the large M/L values natural for old stel- 
lar populations, e.g., Queloz et al. (1995) and Romanowsky 
(2003), who find M/L ratios for old stellar populations as 
high as 7 or 8, our new results arc now at the limit of com- 
patibility with M/L ratios for old naked stellar populations 
for all the well-studied local dSph galaxies. This holds com- 
fortably for all galaxies except Draco, where compatibility 
is found only just at the la level for our lower limits on 
the confidence intervals. Our model thus yields equilibrium 
isothermal configurations able to explain the observed dy- 
namics of the local dSph galaxies, without the need for dark 
matter. 

A consistency check of the above interpretation is avail- 
able from the comparison of the M/L ratios we found and 
the age of the youngest stellar population found in each of 
the galaxies, as inferred from the direct studies of the ob- 
served CMDs of the galaxies in question, given in the last 
column of Table (1). It is reassuring that the highest values 



for our inferred M/L ratios, for Ursa Minor and Draco co- 
incide with the galaxies showing no star formation over the 
past 10 Gyr, the oldest galaxies in the sample, while the 
lowest values for M/L arc obtained for Leo I and Fornax, 
the youngest galaxies in the sample with stellar populations 
as young as 2 Gyr. For the remainder of the sample, galaxies 
showing their youngest stars at intermediate ages, we find 
the intermediate values for M/L. The interpretation of the 
dynamics under the proposed model thus accords with the 
natural increase in the M/L ratios of stellar populations 
due to the ageing of stars and the consequent build up of 
black holes, neutron stars, and white dwarfs. This corres- 
pondence is natural in any modified gravity scheme where 
the stars alone are responsible for the dynamics, but has 
to be thought of as a fortuitous coincidence under the dark 
matter hypothesis. 

As a second consistency check, we now show the res- 
ulting projected surface density mass profile for a model 
of the Leo II galaxy, normalised by the total luminosity of 
that system. We construct the surface brightness profile us- 
ing our model for an equilibrium isothermal system having 
velocity dispersion and projected i?hm (assumed equal to 
the observed i?hi), equal to those observed for Leo II. This 
is given in Figure (2), where we have also plotted the ob- 
servation for the actual surface density light profile of Leo 
II, from the star counts analysis of Coleman et al. (2007), 
out to the radius where measurements fall below the back- 
ground noise level, both normalised to the same total lu- 
minosity. We note that the error bars of Coleman et al. 
(2007) are only a lower estimate of the confidence intervals 
for this comparison, as they only refer to the errors in the 
measured star counts and not to hill surface density profile 
inferences. A very good agreement is evident, showing the 
proposed models are a good, fully self-consistent represent- 
ation of the dynamics and both the integral and spatially 
resolved light distribution in the well-studied local dSph 
galaxies. 



4. Scaling relations 

We now turn to the scalings shown by the dSph galaxies 
in our sample. Firstly, we show in Figure (3) the beha- 
viour of the equilibrium isothermal configurations we are 
solving for, in terms of the resulting projected i?hm as a 
function of the input value of po, at a constant value of 
a = lOkm/s. The remarkable feature of Figure (3) is that, 
after decreasing as po increases, when po reaches a value of 
about 0.3Mopc^'', the resulting _Rhm stops changing and 
settles at relatively constant value of around 150pc, below 
which it does not fall further. This is qualitatively repro- 
duced at all values of a, with only small changes in the 
minimum values of i?hm, for a in the range of values ob- 
served for local dSph galaxies. This is interesting, as it of- 
fers a natural explanation for local dSph galaxies showing a 
minimum lower value for their projected half-light radii of 
around 150pc, with most lying around a factor of 2 above 
this critical limit, as noticed by Gilmore et al. (2007) and 
seen from Table (1). Indeed, all our models constrained by 
the observational values of a and i?hi for the local dSphs, 
occur within the flat region of the -Rhmvs.po space. 

We can now try to understand the assigned Newtonian 
values of M/L for local dSph galaxies and the scalings they 
show with total luminosity or absolute magnitude, e.g.. 
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Figure 3. Dependence of the resulting projected lialf-mass 

radius in kpc, against the assumed central density in 
Mqpc^^, for isothermal equilibrium configurations under 
the proposed MOND formalism, at fixed a = lOkm/s. A 
very broad region where -Rhm remains almost constant is 
evident, at precisely the level found to define the minimum 
iihi values for observed dSph galaxies, which have typical 
values of (7 ~ lOkm/s. 



Mateo (1998). The assigned Newtonian values of M/L will 
never be far from (e.g. Gilmore et al. 2007) 



(M/L) 



N 



hi 



G 



) ■ (10) 



Given our results in the previous section, or alternat- 
ively taking the observed Rhi of around 300pc as an empir- 
ical fact, we can evaluate the above equation to first order 
at i?hi = 0.3. Also, for an average old stellar population, 
we can take a constant M/L value of 5 as a representative 
value for the local dSphs, as our inferences yield, and re- 
place cr^ in equation (9) for the corresponding value through 
equation (10) to yield 



(M/L) 



N 



'375^ 
3.36 



Introducing the absolute magnitude My ~ — 2.51og(Ltot) + 
4.83 and taking the logarithm of the above equation gives 



log(M/L)N = 3.77 + 0.2M^ 



(12) 



Figure (4) now gives a plot of equation (12), superim- 
posed on recent determinations of the Newtonian M/L val- 
ues for a larger sample of local dSph galaxies, using repor- 
ted values of L to calculate the values plotted on the x-axis, 
and of L, a, and Rhi to calculate Newtonian M/L values 
through equation (10). Error bars give extreme la confid- 
ence intervals, including uncertainties in all the parameters 
used. Data as summarised in Walker et al. (2009a) (using 
the revised version of the data for the classical dSphs, as 
given in Walker et. al 2009d), and references therein. The 
good match is evident, extending for over 13 orders of mag- 
nitude in Mv, which is even more remarkable because the 
trend found extends towards the smallest newly discovered 
dSph systems. 



Figure 4. Logarithms of the Newtonian M/L values of 

local dSph galaxies against total observed V band mag- 
nitudes, from data as summarised in Walker et al. (2009a). 
The straight line is not a fit, but gives the Newtonian M/L 
values that would be assigned to isothermal equilibrium 
populations of stars having a constant intrinsic stellar M/ L 
ratio of 5 and a constant half-mass radius of 300kpc, con- 
sistent with what is observed, and to what the proposed 
model yields for equilibrium configurations having velocity 
dispersions of lOkm/s, as the observed dSph show, accord- 
ing to the model presented here. Galaxies in ascending or- 
der of luminosity, with original references: 1- Sgr(23, 24), 
2- Fornax (1,2), 3- And II (17, 18), 4- Leo I (1, 5), 5- Cetus 
(17, 22), 6- Sculptor (1, 2), 7- And XV (20, 21), 8- Leo II 
(1, 6), 9- Tucana (25, 26), 10- Sextans (1, 2), 11- Draco (3, 
4), 12- Carina (1, 2), 13- Cven I (3, 10), 14- Umi (1), 15- 
And IX (19), 16- Leo T (3, 10, 14), 17- Hercules (3, 11), 18- 
Bootes 1 (3, 8), 19- UMa I (3, 8), 20- Leo IV (3, 10), 21- 
Cven II (3, 10), 22- Leo V (12, 13), 23- UMa II (3, 10), 24- 
Coma (3, 10), 25- Willman 1 (3, 8), 26- Bootes 2 (3, 9), 27- 
Segue 2 (16), 28- Segue 1 (3, 15). References: 1- Irwin & 
Hatzidimitriou (1995), 2- Walker et al. (2009c), 3- Martin 
et al. (2008), 4- Walker et al. (2007), 5- Mateo et al. (2008), 
6- Koch et al. (2007b), 7- Martin et al. (2007), 8- Koch et 
al. (2009), 9- Simon & Geha (2007), 10- Aden et al. (2009), 
11- Belokurov et al. (2008), 12- Walker et al. (2009b), 13- 
Irwin et al. (2007), 14- Geha et al. (2009), 15- Belokurov et 
al. (2009), 16- McConnachie & Irwin (2006), 17- Cote et al. 
(1999), 18- Chapman et al. (2005), 19- Ibata et al. (2007), 
20- Letarte et al. (2009), 21- Lewis et al. (2007), 22- Ibata 
& Irwin (1997), 23- Majewski et al. (2003), 24- Saviane et 
al. (1996), 25- Praternah et al. (2009). 



We see that the expectations of the model for the in- 
ferred Newtonian M/L values of local dSph galaxies, which 
are modelled as equilibrium isothermal solutions to equa- 
tion (8) and having typical radii of close to 0.3kpc, agree 
very well with independent measurements. The small scat- 
ter in Figure (4) beyond observational error bars is compat- 
ible with the variations in the actual intrinsic M/L values 
for the individual galaxies, and variations of i?hi around the 
average 0.3kpc used in equation (12), which is not a fit to 
the data, but a first-order estimate within our prescription. 
A very simple explanation for the {M/L) jv values assigned 
to all well-studied dSphs is thus naturally afforded by our 



8 



X. Hernandez^'^ et al.: Dwarf spheroidals in MOND 



assumptions for the most average intrinsic M/L and i?hi 
values, without the need of invoking complex astrophys- 
ical processes (e.g. tidal disruption, strong outflows, etc.) 

or fine-tuning any parameters. To first order, one can un- 
derstand the well-established scalings of figure (4) as the 
extension of the TuUy-Fisher relation, which appears in the 
MOND formalism, down to the dSph regime. 

5. Conclusions 

We used a MOND interpolation function that can be seen 

as the addition of both the Newtonian acceleration and the 
low acceleration limit of the MOND formalism on all scales 
(Bekenstein 2004). For the local dSph galaxies, where the 
application of the MOND formalism has been most contro- 
versial, we show that isothermal cqiulibrium configurations 
characterised by well-defined finite total masses and half- 
mass radii result, giving M/L values in agreement with na- 
ked stellar populations. As a result, no dark matter is now 
needed. The observed scalings in i?hi and assigned [M /L)m 
values as a function of total magnitudes arc explained nat- 
urally, including the prc;viously imrccognised correlation of 
resulting M/L ratios with the relative youth of the stellar 
populations of the individual dSph galaxies, for the sample 
studied. This last is a natural consequence of using the 
MOND prescription, but an odd coincidence under the dark 
matter hypothesis. 

Fixing the optimal ^,{x) for the MOND prescription in 
the dSph regime through the arguments presented here con- 
stitutes an interesting condition on this function in the re- 
gime studied, and it provides a constraint on any global 
IJ.{x) MOND function. 
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